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Abstract 

Multifragmentation in Au+Au collisions is investigated at incident energies 
in the range 100-400 MeV per nucleon by means of a recently developed quantal 
Langevin model. The inclusion of quantum fluctuations enhances the average 
multiplicity of intermediate mass fragments, especially in central collisions. This 
is mainly because the excitation energies of fragments are reduced due to the 
quantal behavior of intrinsic specific heat. 
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1 Introduction 



One of the main goals of heavy-ion physics is to explore the various phases of nuclear 
matter and achieve an understanding of the associated phase transitions. At interme- 
diate energies, the nuclear liquid-gas phase transition is of primary interest, and the 
efforts to understand this phenomenon have intensified following the recent extraction 
of a caloric curve of hot nuclear matter |TJ suggesting that the phase transition is of 
first-order and occurs at a temperature of 4—5 MeV. The first-order nature of the phase 
transition may lead to an enhancement of fragment formation || [| |[, since the rapid 
expansion following an energetic nucleus-nucleus collision causes a transformation from 
a hot to a supercooled gas phase. The resulting dilute system is mechanically unstable 
and fragments are therefore likely to be formed, in a manner similar to the formation 
of water drops. 

The experimental exploration of nuclear phase transitions must rely on nuclear 
collisions, and it is therefore necessary to invoke dynamical models in the interpretation 
of the data. In addition to the complications arising from the non-equilibrium character 
of the collision dynamics, it is also important to ensure that the models employed take 
due account of the quantal nature of the nuclear system. 

In order to illustrate this important point, we note that the nuclear liquid-gas 
phase transition differs significantly from the usual liquid-gas phase transition in macro- 
scopic matter. The largest difference lies in the role played by quantal statistics. In 
the case of usual macroscopic matter, the total energies are to a good approximation 
linear functions of the temperature in both the liquid and gas phases. Thus the effec- 
tive number of degrees of freedom is essentially constant in each phase. In contrast to 
this familiar situation, the liquid phase of a nucleus exhibits an increase in the num- 
ber of activated degrees of freedom as the temperature is raised. For example, the 
excitation energy of a nucleus at low temperature increases like E* = aT 2 (where the 
level density parameter is a w A/(8MeV)), which is a typical quantal behavior, while 
the gas phase is characterized by the usual classical relation E* / A = 3T/2. The two 
curves intersect at T « 12 MeV, which is much higher than the suggested transition 
temperature. This indicates that the quantal statistical nature of the nuclear system 
plays an important role for the phase transition and, presumably, for the associated 
nuclear multifragmentation processes. 

To take approximate account of the quantal features of the evolving nuclear sys- 
tem, we augment the socalled Quantum Molecular Dynamics (QMD) model || @, 0] 
by a stochastic term given by the recently developed quantal Langevin model @, |§. 
The QMD model describes the nuclear many-body system in terms of a product wave- 
function of gaussian wavepackets for the individual nucleons and has been successful 
in accounting for a variety of collision observables, such as collective flow and particle 
production ||. The QMD model and its antisymmetrized versions, Antisymmetrized 
Molecular Dynamics (AMD) |TIJ and Fermionic Molecular Dynamics (FMD) JTTJ, have 
also met with some success in describing fragment production |5|, [7|, [10], [12], [13], |i~4]| . The 
fragments formed in these molecular dynamics calculations are generally far too excited 
to remain intact, so the subsequent decay must be accounted for in order to extract 
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the yield of bound fragments that are measured experimentally. This combination of 
molecular dynamics and subsequent (usually statistical) decay calculation may give 
us a correct picture of fragment formation when the source of fragments is limited 
to one compound nucleus or excited projectile-like and target-like fragments. How- 
ever, when the copious intermediate-mass fragments (IMF) appear simultaneously, the 
above fragment formation mechanism becomes dubious. For example, although the 
calculated IMF multiplicity without the subsequent decay seems to reproduce the data 
reasonably [jT^, [H, the decay processes largely eliminates the IMF yield . 

We suggest that part of the reason for this consistent shortcoming of the entire 
class of QMD-like models may be found in the fact that the model is effectively classical, 
since the equations of motion for the wavepacket parameters have been derived by 
a standard variational principle. The resulting description is then the dynamics of 
the wavepacket centroids and the inherent quantum fluctuations are neglected. The 
presence of quantum fluctuations is signaled by the fact that the wavepackets are not 
energy eigenstates. This basic feature generally causes the expectation value of the 
Hamiltonian to fluctuate in the course of time. 

In order to remedy the situation, a quantal Langevin model was recently devel- 
oped [§, ||. The Langevin force enables the wavepacket system to explore its entire 
energy spectral distribution, rather than being restricted to its average value, and it 
leads to a much improved description of the quantum statistical features. In particular, 
the resulting specific heat now exhibits the characteristic evolution from a quantum 
fluid towards a classical g function of temperature || |J, in contrast to the 

behavior emerging with the usual treatment. Since a change of a fragment's specific 
heat is associated with a change in its statistical weight, the effect is clearly relevant 
for the fragment production problem. 

The key features of the results obtained with the quantal Langevin model are 
the occurrence of larger fluctuations and an enhancement of stable configurations, such 
as bound fragments, as a result of the need to take account of the spectral distortion 
of the wavepackets. The former feature arises from the fact that the wave packets 
are populated according to the strength of the eigencomponent for the given energy, 
and therefore the wavepacket parameters can have larger fluctuations than when the 
expectation value of the energy is conserved. On the other hand, in order to project 
out the appropriate energy component from the wavepacket, it is necessary to take 
account of its internal distortion. The combination of these two basic features is then 
expected to enhance the IMF yield at the final stage: while the larger fluctuations 
allow the system to explore more configurations and thus enhances the yield of primary 
fragments, the latter stabilizes the fragments, since the compensation for the quantum 
distortion effectively acts as a cooling mechanism. 

In our previous work |§, we have applied the quantal Langevin model to a simple 
soluble example and have shown that it leads to the correct quantal microcanonical 
equilibrium. However, in order to apply this model to actual nuclear collisions, the 
treatment of the center-of-mass cluster motion must be suitably modified. In the QMD 
and AMD treatments, the single-particle wave function is represented by a gaussian 
wave packet with a fixed width. Then the center-of-mass wave function of each cluster 



2 



is also a gaussian wave packet and, accordingly, it has zero-point kinetic energy as 
well as an energy fluctuation dependent on the velocity Within the framework of 
the quantal Langevin model, this energy fluctuation causes a quantal fluctuation of 
the CM motion of the cluster even if it is well separated from other clusters and 
nucleons. This phenomenon reflects the expansion of a localized free wave packet and 
should be eliminated, as is usually done in Hartree-Fock and AMD calculations, for 
example, before the binding energy is calculated. In order not to modify the Q-values 
in the fragmentation process, the cluster CM motion must also be removed in the 



dynamical evolution, as was already discussed in refs. fITJ, Therefore, we have 
removed the energy fluctuation of the cluster center-of-mass motions as well. This 
can be accomplished by considering the motion of the individual nucleon constituents 
relative to to the local collective velocity. 

In this letter, we have included the quantal Langevin force into the framework 
of QMD. Although we have formulated the quantal Langevin model within the AMD 
framework, the appearance of the quantal fluctuation originates in the wave-packet 
nature of the many-body state and is therefore a general phenomenon. Thus the basic 
features of the model also apply to QMD which is easier to apply to heavier systems 
and we focus our attention on the fragmentation process in 197 Au+ 197 Au collisions. 
In this system, the total number of nucleons is so large that the statistical properties 
are expected to play a major role. In addition, the multiplicities of intermediate mass 
fragments are measured at incident energies of 100, 250, and 400 MeV per nucleon and 
their impact parameter dependence is also deduced Q. The IMF multiplicities reach 
up to around ten, which is the largest number observed so far. 



2 Quantal Langevin model at fixed temperature 

We first give a condensed description of the recently introduced quantal Langevin model 
for the situation when the system can be regarded as being in thermal equilibrium at 
a given temperature. 

The treatment seeks to take account of the energy fluctuations present in a system 
being described in terms of many-body wave packets. As we have already discussed in 
detail in Ref. |J, this inherent energy dispersion modifies the statistical weight relative 
to the naive the classical form, 

Wp(Z) = (Z\exp(-f3H)\Z) ^ exp(-/3W). (1) 

Here 7i = (Z\H\Z) is the expectation value of the Hamiltonian in the given wave- 
packet state \Z) and thus the last quantity represent the usual classical statistical 
weight. The complex parameter set Z = {z%, z 2 , . . . za} is related to the phase space 
coordinates in QMD, z n = r n /2Ar + ip n /2Ap, where Ar and Ap are widths of wave 
packet. By invoking the harmonic approximation, it is possible to obtain a good 
description of the statistical weight by means of a simple "free energy" , 

T P {Z) = -\ogW p {Z) = % (1 - exp(-/3D)) , (2) 
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where D = o 2 E jE* is the effective level spacing. (The energy of the wave packet relative 
to its ground state is denoted by E* and cr| is the associated variance.) 

The relaxation towards this approximate quantal equilibrium can be described 
by the following Fokker-Planck equation for the distribution <f)(Z) of wave-packet pa- 
rameters, 



D(p 
~Dt 
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(3) 



v< = -E^- (4) 

where represents either r n or p n . It is easy to check that the statistical equilibrium 
distribution, eq = exp(— J-p), is a stationary solution to the above Fokker-Planck 
equation. Moreover, when the classical statistical weight is employed {i.e. when = 
(37i), the drift and the diffusion coefficients of the Fokker-Planck equation satisfy the 
usual Einstein relation, 

Vf = -^M^. (5) 

On the other hand, when the quantal statistical weight obtained with the harmonic 
approximation is used, eq. (0), the relation is modified. For example, if the effective 
level spacing D does not depend strongly on the wave-packet parameters, the drift 
coefficient reduced by the factor a, 

dn l-exp(-/3Z>) 
V^-a^M,,-, „= jLJU.. (6) 

Since a is smaller than unity, the resultant Fokker-Planck equation gives smaller fric- 
tion, thus in effect relatively larger fluctuations will arise. 

It is convenient to solve the Fokker-Planck transport equation by means of a 
Langevin method. Within the framework of QMD the Langevin equation becomes 

p = f-aPMP-v + gP-e, (7) 
r = v + a(3M r ■ f + g r ■ £ r , (8) 

— / = -— (9) 
dp ' dr ' 

M p = g p g p , M r = g r g r . (10) 

Here r and p are the phase-space centroid parameters for the wave packet and £ is 
used to denote random numbers drawn from a normal distribution with two a variance 
equal to two. In these equations, we have omitted the diffusion-induced drift term and 
that part of the mobility tensor that connects r and p. 
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2.1 Energy Fluctuation 



The key quantity in the quantal Langevin treatment is the energy dispersion a\ within 
each wave packet. In Refs.||, [|, the following approximate form was derived, a\. 

° l=& ^- C "-% = ?I Ar2 /n " /n + Ap 2 «„ ■ «„] • (11) 

Here the matrix C arises from the antisymmetrization and, accordingly, it is the unit 
matrix in QMD. There are two problems in the application of the above expression to 
dynamical processes within the framework of QMD, such as nucleus- nucleus collisions. 
The first concerns the spurious zero-point center-of-mass motion of clusters and the 
second arises from the Pauli blocking, as we shall now discuss. 

The energy fluctuation becomes larger when the wave packet moves faster, as 
should be the case, and it results from the zero-point kinetic energy of the fragment 
center of mass. The expression is adequate when the configuration volume is compact 
and all the degrees of freedom can be treated as wave packets. However, as is already 



discussed in Refs. [15fl , this zero-point CM kinetic energy is spurious and should 
be removed from the Hamiltonian. It must then also be removed from the energy 
fluctuation. In addition, it may be unreasonable to carry out the integral over all the 
parameter space near the given value, since no Pauli blocking effect is taken account 
of in the phase space parameters themselves in QMD. Then, it appears more physical 
to reduce the energy fluctuation according to the appropriate Pauli blocking factor. 
In order to satisfy these two requirements we have adopted the following form of the 
energy fluctuation, 



<?l = i 1 ~ /«) 2 Ar 2 f n ■ f n + Ap 2 (v n - ti„) ■ (v n ~ u n ) 



'12) 
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U i = Jq.Y. C H V i ■ ( 13 ) 

1 j 

The additional factor (1 — f n ) 2 represents the Pauli blocking for small changes, where f n 
is the Wigner function. Here Vj denotes the velocity of the nucleon j and the local flow 
velocity it; is obtained by performing a suitable smoothing, using c^- = exp(— r^/4Ar 2 ) 
with Ni = J2j The quantity entering in the energy fluctuation is then v — u, the 
velocity of a wave packet relative to the local collective flow, which represents those 
degrees of freedom having a wave packet nature. 

On the other hand, the collective velocity u, which is thought of as the cluster 
velocity, represents the degrees of freedom having a classical nature. The drift term of 
eq. (0) is therefore modified accordingly, 

p = f - a /3M p -(v-u) - (3M P u + g p ■ £ p . (14) 

It should be noted that when the energy fluctuation a\ is small, the above Langevin 
equation leads to the normal (classical) Langevin equation, p = f — (3M P ■ v + g p ■ £ p . 
This situation is realized when all the fragments are well separated and the intrinsic 
energy fluctuation is very small (i.e. all the fragments are close to their ground state). 
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2.2 Thermal distortion and observation 



In addition to modifying the statistical weight, the energy fluctuation also modifies the 
meaning of wave packet ensemble, since it causes a thermal distortion of the spectral 
strength distribution of the energy eigencomponents within each wave packet. It is 
important to take account of this effect. 

Since the expectation value of an observable O is generally given by 

Tr (6 exp (-/?#)) 

-4 o y = — ^ — . ! J . (15) 

Tr (exp (-/3H)) 

In order to obtain this quantity as an weighted average over the thermal ensemble 
of wave packets, with the weight factor being the usual statistical weight Wp(Z), 
it is necessary to distort the individual wave packets appropriately. The resulting 
contribution to -< O y from a given wave packet is then temperature dependent, 

Op{Z) = (Z\ exp(-pH/2) 6 exp(-pH/2)\Z) . (16) 

The distortion operator exp(— (3H/2) reduces the expectation value of the Hamil- 
tonian in the particular state | Z) . The thermal distortion is calculated by replacing the 
time t by the imaginary time it in the equation of motion. The resulting "evolution" 
is then described by a cooling equation, 

dp n 2Ap 2 dr n 2Ar 2 

-i- = t (v n — u n ) , — = ——/„, (17) 

dr n dr n 

with which the state should be propagated until r = h/3/2. Here, v is again replaced 
by v — u in order to leave the collective (or cluster) motion unaffected. 



2.3 Fragmentation at fixed temperature 



We have incorporated the above described quantal Langevin force into the QMD model. 



In the present implementation, the Gogny interaction [16] is used as the effective inter- 
action. In addition, a correction for the fragment zero-point kinetic energy is included 
as described in refs. [ 10| , [15]. Moreover, since the quantal Langevin force is sensitive 



to the minimum state of the given Hamiltonian, a Pauli potential |15, [17] to guarantee 
that the nuclear sizes and bindings are reasonably reproduced by the model. 

Figure [l] shows the resulting mass distribution of fragments as obtained at various 
temperatures having a low average density (to ensure that the fragments are well sepa- 
rated). The results of the quantal Langevin model agree well with the grandcanonical 
calculation when T > 4 MeV. At lower temperatures, probably because of the limited 
functional space of QMD, both the quantal Langevin model and the classical Langevin 
model fail to describe the population of the heaviest fragments which in turn may lead 
to the overestimation of the proton multiplicity. 
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Nuclear Mass Spectra in Box (p=0.012) 




10 20 30 40 

Fragment Mass 



Figure 1: Mass distribution in a canonical ensemble of 40 nucleons. 

Circles and triangles show the calculated results of QMD at the indicated tem- 
peratures T with (a < 1) and without (a = 1) the quantal modification. As an 
instructive reference for the quantal calculation, we also show the grand canonical 
results (solid lines) based on a simple statistical treatment [18]. 
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It is interesting to note that the mass distribution changes its appearance in the 
temperature region T = 4 ~ 6 MeV. At lower temperatures, a large part of the nucleons 
are bound in large fragments, while at higher temperatures, the lighter fragments 
grow abundant. This behavior reflects the transition from the excitation of various 
collective modes of heavy nuclei to the realization of fragmentation degrees of freedom. 
In our previous work || [|, we have shown that this feature can be also seen from the 
relationship between temperature and excitation energy. The maximum specific heat 
is obtained for the temperature T ~ 4 MeV, primarily because the fragment degrees 
of freedom are activated rapidly as T is increased above this value. 



3 Application to nuclear collisions 



In order to apply the treatment described above to actual nucleus-nucleus collision 
processes, some important modifications of the QMD model are required. In the usual 
QMD treatment, the expectation value of the Hamiltonian is a constant of motion 
(since the equations of motion for the wavepacket parameters have been derived from a 
standard variational principle). When the quantum fluctuations are taken into account, 
the inherent energy fluctuation within each many-body wavepacket makes it necessary 
to include wavepackets having different energy expectation values. At first, this feature 
might seem unreasonable as it appears to violate energy conservation. However, an 
energy eigenstate can be constructed as a superposition of wavepackets whose energy 
expectation values (the diagonal matrix elements of the Hamiltonian operator) are not 
necessarily the same as the energy eigenvalues. This feature grows more prominent 
in the reaction region and there is therefore no reason to exclude wavepackets with 
different energy expectation values in the dynamics, especially during the strongly 
interacting stage. 



3.1 Quantal Fokker- Planck equation at a fixed energy 

By proceeding as in the case considered above where the temperature is given, it is 
possible to derive a quantal Fokker-Planck equation for situations in which the energy 
is specified, such as the evolution of an isolated system. The equilibrium distribution is 
now a microcanonical ensemble, with a corresponding statistical weight for each wave 
packet, 

p E (Z) = (Z\5(E-H)\Z) oc f^^-fM-n/D) . (18) 

The last relation assumes that the spectral function is given by a continuous Poisson 
distribution. The drift coefficient can then be readily calculated, 

Vi = -EM^,-^|, (19) 
F E {Z) = - log p E (Z) , (20) 
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dT E n-E 

l j n = -T^j - 2 ■ [21) 



dH a 



E 



It is interesting to note that the drift term acts to restore the energy: when the Hamil- 
tonian 7i is greater than the given energy E, the drift term reduces the Hamiltonian 
(since fin is positive in this case and the mobility tensor is always positive definite. 
Furthermore, since the phase volume with larger TC is larger, the dynamical trajec- 
tory usually goes through the region with TC > E. Therefore, the drift term generally 
reduces the excitation energy of fragments in the final stage of collision. In the asymp- 
totic region where all the fragments are well separated and all of them are sufficiently 
cold, the energy fluctuation becomes small and the energy expectation value converges 
to the given value. 

Thus, the basic features of the quantal Langevin model remain the same also 
in the case of a fixed energy, namely larger fluctuations and intrinsic distortion. To 
show this important point, we introduce the state-dependent microcanonical inverse 
temperature fl E , whose average over the entire phase space with the weight p E (Z) gives 
the global microcanonical temperature, 

r\ Q \ 

Pe(Z) = — \ogp E {Z) = -— T E {Z) ~ - \og{H/E) . (22) 

The last relation holds when we employ the continuous Poisson distribution and the 
effective level spacing D is small compared to TC and E. In terms of this /3 E , we have 
fin = a E j3 E , with a B — (1 — exp(—j3 E D))//3 E D. Thus the drift term is reduced by 
the factor a E compared to the normal Einstein relation, in analogy with the canonical 
case of a fixed temperature. 



3.2 Mobility tensor 

Although the preceding statistical discussion gives the form of Fokker-Planck equation, 
including the modified fluctuation-dissipation relation, the mobility tensor must be 
determined separately. Although the mobility tensor does not affect the equilibrium 
properties, it is important in a dynamical scenario since it determines the time scale. 
As we discussed in refs. |§ [|, this mobility tensor essentially represents the transition 
rate from the current dynamical state. This transition rate should be described by 
the residual interaction between nucleons, i.e. the off-diagonal part of the original 
Hamiltonian operator. In the present study, we have used the following relatively 
simple mobility tensor, 

T 

M r nhmj = <?<W<% Ar 2 ^ (1 - f n ) , (23) 

M^ mj = g5 nm 5 l3 Ap 2 ^(l-f n ), (24) 

K) 2 = Ar 2 f n -f n , (25) 
«) 2 = Ap 2 (v n - u n ) ■ K - u n ) . (26) 

Since the transition induced by the quantal Langevin force is of one-body nature, we 
have incorporated the blocking factor 1 — f n . 
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Figure 2: IMF multiplicity before statistical decay. 
Circles and triangles indicate QMD results at given energies with and without the 

quantal Langevin force, respectively. The experimental data Q are shown by solid 

diamonds. Dotted lines show the experimental data using a scaled impact parameter 

assuming a maximum impact parameter of 10 fm. The detector efficiency is not taken 

into account in the calculation. 



3.3 Results for Au+Au collisions 

We have applied the quantal Langevin model to collisions of two gold nuclei at inci- 
dent energies of 100, 250, and 400 MeV per nucleon. In addition to the Hamiltonian 
discussed above, an energy-conserving two-body collision term is also included. The av- 
erage multiplicity of intermediate-mass fragments calculated for 197 Au+ 197 Au is shown 
in fig. 0. The treatments with and without the quantal Langevin force agree with each 
other qualitatively. For example, the IMF multiplicity has the peak at central collisions 
at E/A — 100 MeV, while the peak moves to finite impact parameter at higher incident 
energies. 
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Figure 3: IMF multiplicity after statistical decay. 
The display is similar to fig. ^. Circles and triangles show the results of the combined 

QMD and statistical decay model (SDM) calculation with and without the 

quantal Langevin force, respectively. 



Since the fragments extracted from the dynamical calculations are still excited, 
it is necessary to include the subsequent statistical decay before contact with exper- 
iment can be made. Figure [5] shows the resulting average IMF multiplicity after the 
statistical decay chains have deexcited the fragments to below the light-fragment sepa- 
ration energy (we shall consider such relatively cold fragments as being in their ground 
states). As the statistical decay model (SDM), we have used the evaporation model 
of Piihlhofer ]nj. In the case of QMD without quantal Langevin force, the excitation 
energies of the primary fragments are large enough to largely eliminate the IMF compo- 
nent from the final mass distribution. This tendency is more clear for central collisions, 
where the fragments are likely to be highly excited. On the other hand, QMD with 
the quantal Langevin force reduces the excitation energies of the primary fragments 
because of the quantum-statistical nature of the intrinsic degrees of freedom. As a 
result, the difference between the average multiplicities before and after the statistical 
decay is less than one in the quantal Langevin model. 
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Figure 4: Formation time of bound IMFs. 
Circles and triangles indicate the time dependence of the mean number of bound 

intermediate-mass fragments, as a result of the combined QMD and statistical de- 
cay model (SDM) [|l9| calculation with and without the quantal Langevin force, 
respectively, for central collisions (b < 4 fm). 



Although the difference between the two models seems small at E/A — 100 MeV, 
the IMF formation mechanisms are quite different. This feature becomes more clear 
when the ground state population is examined. With the quantal Langevin force, the 
primary fragments emerge relatively cold as a results of the drift term (the intrinsic 
distortion referred to above), while there is no other way than nucleon emission to cool 
fragments in the normal QMD model. Thus, the ground-state IMFs are more easily 
formed in the quantal Langevin model. The evolution of the multiplicity through the 
statistical decay stage is shown in fig. ||. It is clear that a large part of the stable IMFs 
are created in the dynamical stage in the quantal Langevin model, whereas almost all 
the IMFs that appear with the usual QMD treatment owe their stability to the emission 
of particles during in the statistical decay stage. Such a difference can, in principle, be 
subjected to measurement by exploiting the momentum correlation of the fragments. 
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4 Concluding remarks 



In this paper, we have adapted the recently developed quantal Langevin model to sim- 
ulations of nuclear collisions by embedding it into the widely used Quantum Molecular 
Dynamics model in which the individual nucleons are described by gaussian wave pack- 
ets. The essential modification arises from the fact that the nuclear collision proceeds 
in isolation, so it is appropriate to consider the energy rather than the temperature to 
be fixed. 

The essential effect of incorporating the quantal Langevin model is that the sys- 
tem exhibits larger fluctuations and the degree of excitation of the emerging fragments 
is smaller. These two features conspire to enhance the production of intermediate-mass 
fragments. 

The usual molecular-dynamics treatment leads to primary fragments that are 
typically sufficiently excited to emit nucleons. It is therefore essential to add an "after- 
burner" that subjects the unstable fragment to statistical decay. This process causes a 
strong suppression of the IMFs, and a corresponding enhancement of lighter fragments. 
Therefore, even though the usual dynamics leads to a large IMF multiplicity at the end 
of the dynamical stage [|13|, 14j| , the final yields are not necessarily that large ||14|| . On 



the other hand, the quantum Langevin treatment leads to fragments having a relatively 
small degree of excitation and, consequently, a larger proportion of them survive the 
statistical decay. 

The starting point for the present work was the fact that wave packet molecular 
dynamics displays classical statistical features, an observation that has stimulated lively 



discussion |8], |9], [15], [H], |22], |23|. In our previous work |J, we demonstrated that when 
the mean energy is fixed molecular dynamics in a common harmonic oscillator potential 
may give the same observation with and without the quantal Langevin force. However, 
this coincidence depends strongly on the well-bound nature of the system and the 
calculated results in the present paper suggest that additional fluctuations are necessary 
to achieve a satisfactory molecular dynamics description of fragmentation processes. A 
similar conclusion has been reached by Ono and Horiuchi [ 23| . They included additional 



fluctuations to mimic the Vlasov dynamics in AMD, and also performed a constrained 
distortion of the wave function at each point in time. The combined effect of this 
approach seems to be rather similar to the two main features of the present quantal 
Langevin model. 

Before it is possible to quantitatively assess the role played by the quantum 
nature of the system, it is necessary to resolve a number of important ambiguities 
in the quantal Langevin dynamics. For example, the mobility tensor adopted in this 
present work is based more on physical intuition than on first principles. Another 
ambiguity concerns the residual two-body collisions, since the quantal Langevin term 
causes the level density of the final state to differ from the classical one. 

However, taken at face value, our results for Au+Au collisions indicate that the 
inclusion of the quantum fluctuations in the wavepacket dynamics leads to a signifi- 
cant increase in the production of massive fragments with low excitation. Although 
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the experimental data are not yet reproduced quantitatively, the improvement is sig- 
nificant over the results obtained with the usual treatments in which the quantum 
fluctuations are ignored. This general qualitative result suggests that the underlying 
quantum nature of the nuclear many-body system may indeed play a significant role 
in fragmentation reactions. 
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